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Abstract 

We discuss the situation where attractive and repulsive portions 
of the inter-particle potential both contribute significantly to glass 
formation. We introduce the square-well potential as prototypical 
model for this situation, and reject the Baxter as a useful model for 
comparison to experiment on glasses, based on our treatment within 
mode coupling theory. We present explicit result for various well- 
widths, and show that, for narrow wells, there is a useful analytical 
formula that would be suitable for experimentalist working in the field 
of colloidal science. We raise the question as to whether, in a more 
exact treatment, the sticky sphere limit might have an infinite glass 
transition temperature, or a high but finite one. 

KEY WORDS: Colloidal systems, Baxter model. Disordered systems. Glass transition, 
Mode Coupling Theory. 
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1 Introduction 



Kinetic arrest phenomena at a transition temperature Tc occurring just prior 
to the thermodynamic glass transition, Tg, are well known and have been 
studied extensively using the mode-coupling theory and simulation methods. 
Without being exhaustive, we can point to a number of reviews 0] and pa- 
pers 1^, ^, Q in the literature. There has also been considerable experimental 
work on this field [§]-[13|. Typically such phenomena are considered to be 



driven primarily by packing effects where the repulsive part of the potential 
is of primary importance, with the attraction providing merely a modulation 
of the overall phenomena. In this case the hard sphere system exhibits most 
of the relevant phenomena and may be viewed as the prototypical model of 



this type of kinetic arrest [14 



However, there has recently been interest p5| , |T^ in cases where the attractive 
part of the potential, or rather its interplay with repulsion, is more deeply 
implicated in the arrest phenomenon. This can lead to glasses with much 
richer structure, including long-ranged density correlations that are frozen 
into the system. New dynamical phenomena might also be expected in such 
systems, especially where there is subtle interplay between attraction and re- 
pulsion near the glass transition. As in the repulsive case, it is natural to seek 
a prototypical model that can be conveniently studied within mode-coupling 
theory with the aim of elucidating the special features of such a system. It 
is in particular desirable to consider those aspects that represent more than 
simple modulations of hard-core behaviour. One such choice is represented 
by the Baxter model which has been studied by some of the authors 
[p!5| , [T^ in previous publications using the Percus-Yevick approximation 



for the static structure factor, and the mode-coupling theory. The results 
are promising in that novel dynamical phenomena do emerge, and also it 
is possible to create kinetic arrest where a long correlation length scale is 
quenched into the system. 

Of course, some of the limitations of a model with hard-core potential and 
delta-function like attractive potential, as treated by Percus-Yevick ||18|] are 
evident. For example, it is believed that such a potential would lead only to 
disordered and crystalline phases, rather than the liquid-gas phase separa- 
tion implied by the P-Y approximation [T^. On the other hand, it may still 
be possible to make some progress with this approach. For example, in the 
disordered phase the P-Y approximation is reasonably good for the structure 
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factor so in those cases where the kinetic arrest arises, prior in temper- 
ature or time to a crystalhzation, we may expect reasonable results. It may 
even be the case that the capability to study arrest near what is believed 
to be a 'metastable' liquid-gas transition may be of some practical interest 



21|, ^ |2J]. Here we discuss some unsatisfactory features originating 
from the large-g-tail of the static structure factor, in conjunction with the 
study of the ideal glass transition in the frame of the mode coupling theory 
(MCT). 

It is perhaps prudent to mention at this stage that the use of Mode Coupling 
Theory (MCT) for such questions is itself not without controversial aspects. 
There are numerous criticisms of MCT, and in any case, even if one does 
accept it as a viable mean-field-type theory, there remains the criticism that 



it neglects such effects as 'hopping' and other phenomena [^. These may 
be relevant in square-well systems. 

However, we shall further argue that, within the MCT, the P-Y approxima- 
tion and possibly even the exact structure factor of the Baxter model, leads 
to very high or infinite glass transition temperature for a broad range of 
packing fractions in the region where attractive interactions dominate. The 
outcome of our deliberations shall be that we shall propose as our proto- 
typical model the 'narrow' square well potential, treated by P-Y and MCT. 
We shall show that, for this case, the structure factor is damped at large-g 
due to the finite well width, and that this leads to finite values of Tc, the 
effective transition temperature, throughout the phase-diagram. Moreover, 
the 'narrow' square well potential seems to share most of the feature of the 



kinetic arrest discussed in earlier works [|T5|, |16l. We shall conclude that the 
underlying phenomena reported there are therefore robust and likely to be 
observed experimentally, even though the Baxter model is flawed. 
The organization of these observations is as follows. We shall first consider 
the static structure factor of the square well potential as the input to the 
MCT. This will permit us to explain new features introduced by a finite 
well-width. Next we shall consider the effects of these features on the MCT 
prediction of the transition. Finally we shall discuss the general behaviour 
of the arrest phenomenon as a function of well-width. 
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2 The P-Y Approximation to the Square- Well 
Potential 



We consider a square well potential with repulsive core of diameter R and 
total range R' . The well width is therefore parametrized by e = , and 
the well depth by u. Note that by taking the appropriate limit of well depth 
to width we may recover the Baxter model which possesses an analytical 
solution for the structure factor. By the Baxter sticky sphere model we 
mean e — and the well depth becomes infinite according to ■^ln(12re), r 
being an effective temperature and (3 = . To solve the Baxter model 
within PY an external parameter A is introduced to define h{r) within the 
range of the attraction. It is then easily calculated through a second order 
equation arising from the matching conditions that have to be verified. Since 
the limit A — >■ corresponds to the hard sphere limit, it may be thought as 
an attractive energy scale. 

Although no analytic solution exists for the square well problem, this model 
potential is simple enough to be exhaustively studied. Thus the integral equa- 
tion for the square well problem is readily derived using Baxter's approach. 
The equations are 

rc(r) = -Q'{r) + 2np dtQ'{t)Q{t - r) (1) 
for < r < R', and 

rh{r) = -Q\r) + 27ip T dt{r - t){h{\r - t\)Q{t) (2) 

Jo 

for r > 0, where c(r) is the direct correlation function, h{r) is the indirect 
correlation function and Q' is the derivative of Q{r) The function Q{r) 
is related to S{q) via 

Siq)-' = Q{q)Q{-q) (3) 

and 

Q{q) = 1-271 p dre^^'Qir). (4) 
Jo 

This equation can be solved numerically using the P-Y approximation as clo- 
sure, and it is via this expression that the potential enters the formulation. 
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Of course, the PY closure has many hmitations. However, it is quite accept- 
able for hard spheres, and it is expected to be reasonable for short-ranged 
potentials . 

Instead of a numerical solution, the P-Y theory may be solved in a series 
expansion for small well-width, e,a program that was partially carried out by 
Menon et al. ||2^. The result for the structure factor is. 



S{q) 



1 = 24r/ 



2r]X r 



/2(e?) - 2/3(^9) 



+ 



- 24r/[/2(g-)-(l-e)V2((l-e)g); 
where q = qR', = f pR'^ and a, (3 and /i are given by 

(1-277-^)2 



(5) 



a 



3r/(2 + rif - 2^(1 + 77/ + rj") + p^{\ + 77) 



2(1-77)4 



(6) 
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{4x^ — 24x) sinx 



[x 



12x2 + 24) cosx + 24 



(7) 



a formula that was derived by Liu et al. Note that 77 is the effective 

volume fraction calculated according to the large length, R', that defines 
the full range of the potential. Sometimes we shall also use the real volume 
fraction calculated with R, the hard core diameter, and this is named <p. The 
two quantities are the same in the limit of Baxter spheres. 
In fact, one can push this asymptotic theory a little further than has been 
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done previously. That is, we can solve all parts of the theory to order e, 
including the A-equation, obtaining the following, 





- 12A 








+ 



+ 



where the physical solution A(e) is now of the form: 



(8) 




(9) 



It is fairly straightforward then to compare the results of equation where 
A = A(e) with numerical solution of equations(0, H, ^, It is also possi- 
ble to illustrate differences between the Baxter sticky sphere and the square 
well model. We shall here present a representative selection of results from 
three values of e (e = 0.01,e = 0.03 and e = 0.09), chosen to illustrate var- 
ious aspects of the mode-coupling calculations. The smallest value of e is 
representative of sticky- sphere type behaviour, but with a finite temperature 
transition. The largest value (e=0.09) already exhibits properties closer to 
that of a square well. We note first that, for e = 0.01, the agreement between 
the numerical calculation of equations (|1], 0, ^ ^) and the leading order in the 
well- width expansion is essentially perfect (Fig.||), and we may accept equa- 
tion as defining the P-Y approximation to the small well-width problem. 
As one might expect for such small e, the moderately small-g (corresponding 
to the nearest neighbour distances) behaviour of the structure factor in equa- 
tion (^) is in good agreement also with the Baxter model. However there are 
some very important differences for the large-g behaviour of S{q) as illus- 
trated in Fig. ^ There are two points to note here that relate to later MCT 
calculations. Firstly the MCT calculations are sufficiently delicate at large-g 
for it to be necessary to include many momentum points on the numerical 
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Figure 1: a) Square well numerical result plotted versus e-expansion fore = 
0.01. b) Detail of the first peak of the structure factor (the dotted line is the 
numerical calculated square well, the continuous line is the e-expansion). c) 
Detail of the small-q region of the structure factor. 



grid. It is therefore highly advantageous to have a formula such as equation 
@, rather than a numerical solution. Secondly from Fig.Q (and later com- 
parison of Fig.^) we see that beyond ^a, the structure factor of the square 
well model decreases more rapidly than that of the Baxter model. In fact 
consideration of finite well-width leads, for many purposes, to the effective 
introduction of a cut-off at twice the well width. 



It is this reduction of the range of S{q) that leads ultimately to meaning- 
ful mode-coupling calculation. 

Further comparisons between the e-expansion and the numerical results from 
equation (|l|) are also possible. As might be expected, at e = 0.09 the devi- 
ation between the two results becomes significant (Fig.|^), both for small-g 
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Figure 2: Structure factor at large-g of the Baxter model and e-expansion 
for e = 0.01. The value qx=-^ and by the examination of Fig.^ it is evident 
that much beyond this point there is little contribution to the transition 
temperature. 



and the first peak, and thereafter one must rely on the numerical solution. 
Even so, it is interesting to note that the large-g behaviour is still quite well 
represented by equation (J^). 

We now turn to the thermodynamic aspects of the P-Y solution of the finite 
well problem. We have noted earlier some intrinsic failings of the P-Y ap- 
proximation of the Baxter model in that liquid-gas rather than the expected 
crystallization transition is produced [|1^. There are other unusual features 
of the P-Y solution, such as the asymmetry of the phase diagram with respect 



to the spinodal (see for example jS^, ^]) and the uncertain status of the 



phase diagram to the left of the critical point. In fact, a careful analysis of 
the expansion in square- well width leads to a few new features, not present in 
the Baxter solution. Amongst them there is the observation that the phase 
diagram is no longer so asymmetric with respect to the critical point, and 
there is a finite region of real solution beneath the spinodal curve, near the 
critical point [^. These details need not concern us here, since we shall be 
interested in points of the phase diagram above the P-Y phase separation. 
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Figure 3: Square well versus e-expansion for e = 0.09. 

For our present purpose it is possible to use the P-Y structure factors for tem- 
peratures and densities to the stable (right) side of the high density branch 
of the spinodal line. 

In summation, then, we believe that the P-Y approximation, and in partic- 
ularly the leading order of the small well-width expansion given by equation 
(H) is, at least, a useful approximation to the structure factor for an inter- 
esting range of parameters of the square well potential. We shall now study 
thus range within MCT. 




3 The Mode Coupling Theory 

The mode coupling theory provides a description of the kinetic arrest phe- 
nomenon. For certain values of the interaction parameters the density- 
density correlation function 

possesses a long-time decay with a non-zero infinite time limit. When this 
occurs, diffusion slow dramatically and the viscosity diverges: the glass tran- 
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sition occurs. The Zwanzig and Mori ||3^ formalism and MCT ideas |1| lead 
to the equations 



+ Ql(l),{t) + ly^t) + nl mqit - t')(t)q{t')dt' = (12) 
where mq{t) is given by 

^1=1 J ^^(q'k)0fc(t)0|q_k|(t) (13) 

the vertex function is 

y(q, k) = ^ (q • (q - k)c|q_k| + q ■ kctf SgSkS\^^k\ (14) 

and the two quantity Qq and Uq are respectively the characteristic frequency 
and a white noise term due to the fast part of the memory function; they are 
defined as, 

n = ^'^^^ 

mS{q) 
Uq = Viq^ 

and z/i = 1 in our calculations. In the limit t ^ oo we find the equation for 
the static problem 

/, 1 r(Pk 



1-fq 2 J (27r) 
where fg is the Edwards-Anderson factor: 

fq=hmcPq{t), (16) 

Sq is the static structure factor and Cq = -^g- is the Fourier transform of the 
direct correlation function following the Ornstein-Zernike relation. 
It is clear that /g = is always a solution of the equation (|l^). In fact 
it is possible to show [Q that, where the potential is purely repulsive and 
at low densities, this is the stable solution, and the system is in the liquid 
phase. Above a certain (p^ (that depends on the parameters of the system 
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) non-zero solutions begin to appear and non-ergodic behaviour appears. 
The ideal glass line, moreover, can be defined studying the behaviour of the 
eigenvalues of the stability matrix of the system Here we have bracketed 
the transition line by iteratively solving equation (0). Where more precise 
results have been required we have studied the stability matrix. Of course, all 
these calculations are subject to numerical error, particularly with respect to 
the numerical momentum cut-off in the integrals. Since the attractive glass 
requires consideration of both repulsive and attractive parts of the potential, 
the former length scale being set to be the well-width, we have carried out 
finite-size analysis of our results to check their validity. Now we are in a 
position to solve the kinetic arrest temperature, function of well 

parameters and the structure factor. From now on we shall refer to the 
effective temperature r as introduced by Baxter and related (for small e) to 
the temperature via: r = Y^^^^fc^T WJ^' '^^■^^ -^^^^ point to note is that, 
for the Baxter model in the region where attractions are relevant (previously 
labelled as B2 in references |15| and |]l6l), there are, strictly speaking, very 
substantial contributions from the long tail of the Baxter solution, not taken 
into account in earlier works. The apparent transition observed previously 
is a consequence of truncating the long tail of the Baxter solution. Indeed, 
we find that Tc rises with the number of g-vectors taken into account in the 
calculation. 

For example, in the case of Fig.|| (e=0.01 and 0=0.3881) we see that the 
continues to change with the number of points in the numerical grid until we 
reach a well- width cut-off equal to twice the well width (^).Of course, the 
cut-off required is somewhat dependent on the property of the system chosen. 
However the general outline of the observation is clear. The attractive glass 
transition is strongly finite-size dependent due to the contribution to the 
integral equation (plSf), from momentum scales up to around Q\=-^- Thus, 
for very small well widths the calculations require many thousands of points, 
and for the Baxter model there appears to be no useful number of points 
that can ensure convergence of the transition temperature in the attractive 
region of the transition. The implication is clear. It is the (small) finite- 
well-width model that is the appropriate prototype of the system we discuss 
rather then the literal Baxter Model. However, we make another observation, 
illustrated quite well by Fig.|[ That is, paradoxically, the Baxter model 
solution can be applied in the mode coupling theory to study attractive 
glasses providing the appropriate well-width cut-off {q\) is applied. We can 
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Figure 4: Behaviour of Tc in function of the cut-off as evaluated from equation 
(|15D. e=0.01 and 0=0.3881. As before qx = ^. 



see this because Tc from the Baxter solution truncated at qx is comparable to 
Tc from the true P-Y square well potential of width e = In some cases 
in the literature we have therefore the interesting situation that the Baxter 
structure factor appears to work due to numerical cut-off, even though the 
full model has an unrealistic glass transition temperature. The discussion is 
not purely academic since experimentalist fit data to such models. In future 
we recommend that equation (j^) be applied up to around e=0.03 for the case 
of attractive colloidal particles. Beyond this value of e, more care must be 
taken. 

We now turn to another aspect of this question. We have calculated the 
asymptotic function of e for (j)=0.2 (Fig.^. Evidently, as e decreases 

the glass temperature rises sharply (due to the large-g tail) so that at e=0.001 
the effective temperature Tc is very high, roughly seventeen times higher than 
the Baxter gas-liquid critical temperature. For numerical calculation, we are 
unable to conclusively show the asymptotic limit Tc (e — * 0). However, 
we believe that the glass appears to be the stable solution for all effective 
temperatures where the system does not behave as a herd-sphere fluid. This 
is quite a practical question. Colloidal particles with vanishing well-width, 
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Figure 5: Behaviour of t^. in function of e at (^=0.2 



but large well-depths would therefore not undergo the expected fluid-solid 
T9(] , but would generally found to be a glass, up to temperature 



transition 



where attraction are irrelevant in any case, and the particles behave as hard 
spheres. Since, for finite e, we now have finite values of Tc, it is natural to 
ask if the important dynamical features of the model described in reference 



15] are recovered. The existence of a logarithmic decay of density-density 



correlation, associated to the presence of a cusp singularity in the parameter 
space in the case of "narrow" square well potentials is under investigation. 



4 Conclusions 

As we have noted earlier, it is of interest to identify a simple model, and 
treatment of that model that exhibits the principal phenomena where the 
attractive part of the potential is important in glass formation. We have 
seen that one good candidate is the square well potential where, for very 
narrow wells, important new properties are already evident. Thus, glasses 
with frozen density fluctuations, including those where the correlation lengths 
^ is quite large, appear beneath the attractive glass curve. In addition we 
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have seen that well-width is quite an important parameter. Truly sticky 
spheres in the Baxter sense, or experimental approximations to these, prob- 
ably are nearly always found as glasses, perhaps the powdery precipitates 
found in particle glasses being examples of this. For more moderate square- 
well parameters we may expect interesting glass behaviour, and the details 
of the transition and attendant dynamics are under examination. We have 
also noted that, from the experimental point of view, the finite- well structure 
factor formula (^) represents a more realistic approach to the dynamics of 
such system, though indication of this were already clear much earlier |^ . 



Finally, we offer the conjecture that the exact model of Baxter spheres would 
be a glass for any effective temperature which preserves the presence of the 
attractive interactions in the potential. This could be in agreement with 
Stell's clustering ideas for the Baxter model. In addition, however, we 
have illustrated how the situation evolves away from the Baxter limit. Proof, 
or modification, of the Baxter limit scenario may be an interesting challenge 
to some of the readers of this paper. 
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